Whittaker functions and related stochastic processes 
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Abstract. We review some recent results on connections between Brownian motion, Whittaker functions, 
random matrices and representation theory. 

1. Harish-Chandra formula, Duistermaat-Heckman measure and Gelfand-Tsetlin patterns 

Define J\(x) = h(X)^ 1 det(e XiXj ), where h(X) = Yii<j(^i ~~ A;)- F° r each x, J\(x) is an analytic function 
of A; in particular, Jq(x) = (jYjZi 3^ h(x). The functions J\(x) play a central role in random matrix theory. 
For example, if A and X are Hermitian matrices with eigenvalues given by A and x, respectively, then 

(1) / e tTAUXU^ du= Jx(x) 



IU(n) Jo(x) 

where the integral is with respect to normalised Haar measure on the unitary group. This is known as the 
Harish-Chandra, or Itzykson-Zuber, formula. 

Let j3 = (j3t, t > 0) be a standard Brownian motion in M. n with drift A. Denote by V x the law of /3 
started at x and by the corresponding expectation. Set 

Q = {x e M" : xi> x 2 > ■■■> x n }, T = mf{t > : & £ Q}. 

For \,x € R", write \{x) = XiXi. 

Proposition 1. Forx,\en, J x (x) = h(X)- 1 e x ^F x (T = oo). 

Proof. This is well known, see for example [6]. The function u[x) = V X (T = oo), x € fl, satisfies 
|A« + A • Vu = 0, vanishes on the boundary of fl and lim^^oo u(x) — 1. Here we write x — > oo to mean 
Xi — Xi + i oo for i — 1, . . .,n— 1. Hence v(x) = e xlyX ^u(x) satisfies At; = Xfv, vanishes on the boundary 
of fi and lim^^oo e^ x ^v{x) = 1. The function det(e XiXj ) also has these properties, so by uniqueness, 
v(x) — det(e XiXj ), as required. □ 

The Harish-Chandra formula has the following interpretation. Pick U at random according to the 
normalised Haar measure on U(n) and let fj, x (dy) denote the law of the diagonal of the random matrix 
UXU*. Then the integral becomes 

e ixAuxv dU = [ e A <"V(dl/)- 

U(n) JR" 

Setting m x {dy) = J (x)[i x (dy), we obtain 

e x ^m x (dy) = J x {x). 



The measure m x is known as the Duistermaat-Heckman measure associated with the point x S fi. It has 
the following properties, which are well-known. The symmetric group S n acts naturally on K™ by permuting 
coordinates. The support of the measure m x is the convex hull of the set of images of x under the action of 
S n . It has a piecewise polynomial density. This comes from the fact, which we will now explain, that the 
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Duistcrmaat-Heckman measure is the push-forward via an affine map of the Lebesgue measure on a higher 
dimensional polytope known as the Gelfand-Tsetlin polytope. 

Let x € n and denote by GT(x) the polytope of Gelfand-Tsetlin patterns with bottom row equal to x: 

GT(x) = {P ktj , l<j<k<n: P Kj+1 < P k ^ ltj < P Kj , 1 < j < k < n, P n , = x}. 

Define the type of a pattern P to be the vector 

(n n—1 \ 

^1,1,^2,1 + p 2 , 2 - Pi,i, . . . ,£ - E Pn -^ ■ 
3=1 3=1 / 

Consider the map from U{n) to GT(x) defined by U i— > P where: for each 1 < k < n, Pk t . is the vector of 
eigenvalues of the k th principal minor of UXU* . It is well-known (see, for example, [2] or [I] Section 5.6] for 
a more general statement) that the push-forward of Haar measure under this map is the standard Euclidean 
measure on the polytope GT(x). Moreover, the diagonal of the matrix UXU* is equal to the type of the 
pattern P. From this we obtain another integral representation for the function J\ as 

(3) J x (x) = f e xt ^ eP dP. 

JGT{x) 

2. Whittaker functions 

Set H = A — 2 Y^i=i e~ ai ( x \ where on — ei — ej+i, i = 1, . . . ,n— 1. Write H = for the moment; we 
will drop the superscript again later, whenever it is unnecessary. For convenience we define = d 2 /dx 2 
and ^\x) = e Xx . Following [18], for n > 2 and 9 E C, define a kernel on R" x R" _1 by 

(/ n ri-l \ Ti-1 \ 

9 E x * - £ v* - £ ( eV *~ x * + e Xi+i ~ m ) ■ 
\i=l 7=1 / 7=1 / 

(n) 

Denote the corresponding integral operator by Qg , defined on a suitable class of functions by 

Q e n) f(x) = J Q^(x,y)f(y)dy. 

in) 

The Whittaker functions ip x , A £ C n are defined recursively by 



As observed in [18] . the following intertwining relation holds: 

(5) - 9 2 ) o = o H^-V. 

This follows from the identity —9 2 )Q < "J l \x, y) — Hy l ~ 1 'Qg(x, y), which is readily verified. Combining 

(UJ with the intertwining relation ([5]) yields the eigenvalue equation: 

(6) H^ = (p\l\^. 

Let us now drop the superscripts and write H = H^ n \ ?/>a = W\ • Iterating @ gives the following integral 
formula, due to Givental [23], [28] [18] 

„ n — 1 k 

(7) i>x{x)= ^nn^ 

Mx) t=i f=i 

where T(x) denotes the set of real triangular arrays (X/^, 1 < i < k < n) with T n ^ — Xi, 1 < i < n, and 

n / k k—1 \ n — 1 k 



^a(T) = £ A fc £ T M - £ T fc _ M - £ £ (e—— 

fe=l \i=l 7=1 / fc = l 7=1 
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Now, it is shown in [3] that, for each A € f2, the equation H f = A?/ has a unique solution / = fx 
such that e~ x ( x > f\(x) is bounded and \\m x ^. +00 e~ xtyX ^ f\(x) — 1, where we write x — > +00 to mean ai(x) = 
Xi — Xi+\ — > +00 for each i. Moreover, by Feynman-Kac, 



(8) fx(x) = e x ^E x 



exp 



where /3 S is a Brownian motion in R™ with drift A as in the previous section. The relation between the 
functions fx and the Whittaker functions ipx is thus determined by the following proposition. 

Proposition 2. For A e fi, 

(9) Um e- x ^^Jx(x)=T]^(^-^)- 

PROOF. We prove this by induction on n using the recursion (Q}. Write -0a = ip^ as before, setting 
ip£\x) = e Xx . Then e^Vl^) = 1 and, for n > 2, 

„ / n / 71 n— 1 \ n— 1 \ 

e - A ^)^ n) (x) = / exp - V Ai^i + A n V Xi - V W - V (e"«-«' + e"+*-*<) 
•^R- 1 V i=i \i=i i=i / i=i / 

x V , i"~ 1 i ) A„_ 1 (yi, ■ • -,yn-i)dyi . ..dy n -x 

„ / n— 1 n — 1 \ 

= / e Er= 1 1 (A ! -A„)y l exp _V e» - V e^ 1 "*^ 

R" ~ 1 V i=i *=i / 

x s "=i Al(xi+y * Vi"~ 1 i ) A„_ 1 (a;i + Vi, ■ ■ ■ , x n -i + y n -x)dyx . . . dy n - X . 

By induction, we immediately conclude that, for each n, if k, A G O then e~ x ^ x ' ] tp x n '(x) < Yl i< j r(A^ — Xj). 
Here we are using 

« / ri— 1 n—1 \ 

/ eE7^(M-K)yi exp - V e Vi - V e'^-t-v* ) d Vl . . . dy n - X 

„ / n— 1 \ n— 1 

< / i e^i 1 ^-^)" exp I - £ e» . . . = [] T(A 4 - A„). 

R™ \ i— 1 / i— 1 

It follows, again by induction and using the dominated convergence theorem, that (J9)) holds for A e SI. □ 
Corollary 1. For A e O, 

-E / e^^ds . 

Corollary 2. for i,Ae!l, 

Jx(x) = lim /3-"("- 1 )/ 2 ^(^). 

Proof. By Proposition [TJ the statement is equivalent to 

lim (3- n(n - 1)/2 ^x/p(M = h(A)-VWp,(T - oo). 

0— >oc 

This follows directly from (|10[) by Brownian re-scaling. □ 

As shown in [3], the function ipx, which can be defined by (flU)) . is a class-one Whittaker function, as 
defined by Jacquet |27] and Hashizume [25 . In the notation of the paper [3] we are taking n = {cti/2, i — 
1, ,..,n — 1}, m(2a) = 0, |ry Q | 2 = 1 and tj) v (x) — 2 q k u (x) where q — n{n — l)/2. In the paper [21] , the 
relationship between Givental integral formula and a recursive integral formula due to Stade [52 based on 
Jacquet's definition (see also [26]) is described. 
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Givental's integral formula ([7]) has a very similar structure to the formula ([3]). Indeed, if we define the 
type of an array (T^i, 1 < i < k < n) to be the vector 

(n n— 1 

Ji,i,t 2) i + r 2 ,2 - . . . , 2jr n j - 2J T n -ij J . 



and a measure 



1 — 1 k n—lk 



9 (dT) = n n e M e i+i dT ^ = eTo(T) n u dT ^ 

k=l i=l fc=l i=l 

then 

<M*) = / e^^dT). 
■/r(z) 

On the other hand, if we replace the functions e _e " in the reference measure g by the indicator functions 
^x<y to get a new reference measure 

n—1 k 

g {dT) = Y]_ Y\. lT k^<T k+ i,, l T k+1 , i+1 <T k , t , 

k=l i=l 

then ([3]) can be written as 

Jx(x) = / e xt ^ eT g {dT). 
Jt(x) 

We note the following. If A G iR" then ^ A (a;) = V-a(»); if A € iR" and z/ e 1", then \ip x+l/ (x)\ < ip v (x). 
For each x £ R n , ip^i x ) is an entire, symmetric function of A e C n |21L 125), I31j . There is a Plancherel 
theorem |55L I49|, \21\ 131] which states that the integral transform 

(11) /(A) = / f(x)i> x (x)dx 



is an isometry from L2(R™, dx) onto L^^iR™, s„(A)dA), where L^ m is the space of L2 functions which are 
symmetric in their variables, 1 — \f— T and s„(A)c?A is the Sklyanin measure defined by 

(12) s n (X) = — ±— IJ^-Ak)- 1 . 

For x,[i £ R™, denote by tr^ the probability measure on the set of real triangular arrays (Tk t i)i<i<k<n 
defined by 

/ fda; = M^r 1 [ KT)e^ {T) nn^' 

3 M*) k=iti 

Define a probability measure 7^ by 

The probability measure j x — 7g is the analogue of the (normalised) Duistermaat-Heckman measure in 
this setting. The integral operator K with kernel K(x,dy) — ifio(x)j x (dy) satisfies the intertwining relation 
HK = KA. We can write K(x, dy) = k(x, y)p x (dy), where k is a smooth kernel from R n to R™ = {y £ R n : 
J2iUi — J2i x i\ an d Px denotes the Euclidean measure on R™. For n = 2, 

k(x,y) =exp(-e X2 ~ Vl - e Vl ~ Xl ) 

and, for n = 3, 

k(x,y) = ^\a,b) = 2X (2e( b - a )/ 2 ) 

where 

e -a = e x s -yi- ya +e -xi ) e 6 = e 9 

and -Ko denotes the Macdonald function with index 0. 
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3. Interpretation of "f x in terms of Brownian motion 

A reduced decomposition of an element w € S n is a minimal expression of w as a product of adjacent 
transpositions, that is, w — . . . Sj r , where s; denotes the transposition (z, i + 1). We will also refer to 
the word i = i\i2 ■ ■ ■ i r as a reduced decomposition. By definition, any reduced decomposition has the same 
length l(w), defined to be the length of w. There is a unique longest element in S n , namely the permutation 

/ 1 2 ■•• N ' 
Wo= [n N-l ■■■ 1 

Its length is n(n — l)/2, as can be seen by taking the reduced decomposition 

i = 1 21 321 ... nn-1 ...21. 

The symmetric group acts on R n by permutation of coordinates, and as such is an example of a finite 
reflection group. It is generated by the hyperplane reflections = s ai , i = 1, . . . , n — 1, defined for x € K™ 

by 

where on = e, — e,-+i. Note that Sj corresponds to the adjacent transposition (i,i + 1). 
For a continuous path n : (0, oo) — > R™, define Tj = T tti by 



T,7 7 (t) = jtft) + (^log^ e- Q -^ s »d^ a,, * > 0. 



Let w = Si t ■ ■ ■ Si r be a reduced decomposition. Then 

T — T ■ ■ ■ T 

depends only on w, not on the chosen decomposition [6] . 

We now introduce a probability measure P under which 77 is a Brownian motion in M" with a drift /j, 
and 77(0) = 0. In this setting, a very special role is played by the transform T^ n > = T wo . In the following we 
use the fact that this is well-defined for each n. Write n = (rj 1 , . . . , rf 1 ). For each k < n, set 

(T ktl ,...,T k<k )=TW(r 1 1 ,...,r, k ). 

The evolution of the triangular array T kl j, 1 < j < k < n is given recursively as follows: dT\_\ = dr) 1 
and, for k > 2, 

dT M = dT k _ ia + e T ^- Tk ~^dt 

dT ka = dT fe _i. 2 + ( e T >*.3- T ''-u - e r*, a -r fc _i,i\ ^ 



dT fcjfc _i = dT fe _ life _i + ( e T *=.fc- r fc-i.fc-i _ e Tfc.*-i-Tfc-i,fc-2) ^ 
(13) dT fc! fc = dr7 fc - e Tk - k ~ Tk ^ 1 - k - 1 dt. 

The process, which is clearly Markov, contains a number of projections which are also Markov. For example, 
setting = T k ^ k , we have, for k < n, 

d£ k = dn k - e^-^-'dt. 

This defines a simple interacting particle system on the real line, which has very nice properties. For example, 
in the coordinates and — £j, 1 < i < n — 1, it has a product form invariant measure, that is, a 

product measure which is invariant. 

A remarkable fact is that each row in the pattern T k j is a Markov process with respect to its own 
filtration. This gives an interpretation of the measures 7^ and er^ defined in the previous section. 

Theorem 1. 41 T Wo r)(t), t > is a diffusion process in M. n with infinitesimal generator 

^ = \% Y (i?-f]^ 2 )V'M = ^A + Vlog^.V. 




1 I I — I 1 1 1 1— - 

tl t2 tz tN-2tN-l t 

Figure 1. An up/right path cf> = {0 < ti < .. . < t n -i < t}. 

For each t>0, the conditional law of {Tk^(t), 1 < j < k < n}, given {T Wo r](s), s < t; T Wo rj{t) — x}, is a x 
and the conditional law of T)(t), given {T Wo ri(s), s < t; T Wo rj(t) = x}, is 7* The law of T Wo rj(t) is given by 

v£{dx) = e-^^ t/2 tf;^(x)9t(x)dx, 

where 

(14) 9 t (x)= ( ^x(x)e^i x * t / 2 s n (X)dX. 

J M. n 

In the case n = 2, this is equivalent to a theorem of Matsumoto and Yor [3TJ . 

Write C = Cq and v t — . The diffusion with generator L is the analogue of Dyson's Brownian motion in 
this setting and the measures v t and 9 t (the latter requires normalization) are the analogues of the Gaussian 
unitary and Gaussian orthogonal ensembles, respectively. The diffusion with generator £ M was introduced 
in [3). When fj, € CI, it can be interpreted as a Brownian motion in M™ killed according to the potential 
J2i e Xi+1 ~ Xi and then conditioned to survive forever |29L 130] . The path-transformation T W(I is closely related 
to the geometric (lifting of the) RSK correspondence introduced by A.N. Kirillov [34 and studied further by 
Noumi and Yamada [401 . A discrete-time version of the above theorem, which works directly in the setting 
of the geometric RSK correspondence, is given in [14] . In the discrete-time setting the Whittaker functions 
continue to play a central role. See also [HI 1 1 L 1131 1241 1431 146] for further related developments. 

4. Application to random polymers 

The following model was introduced in [44] . The environment is given by a sequence B\ , B2, ■ ■ ■ inde- 
pendent standard 1-dim Brownian motions. For up/right paths <f> = {0 < t\ < . . . < ijv-i < t} (as shown in 
Figure 1), define 

E(cj>) = B^h) + B 2 (t 2 ) - B 2 {t x ) + ■■■ + B N {t) - B N (t N ^), 

P{d<t>) = Z^rV^'c^, Z?(J3) = J e 0E ^d4>. 

Set X r {{t) = log Z? and, for k = 2, . . . , n, 

X?(t) + ---+X r k l {t)=\og J e B ^ + - +E ^d<t>x...d<i>k, 

where the integral is over non-intersecting paths <f>i, . . . , (f>k from (0, 1), . . . , (0, k) to (t, n — k + 1), . . . , (t, n). 
Let r] — (B n , . . . , B{). Then X — T wa -q and the following holds. 

Theorem 2. [41] The process X(t),t > is a diffusion in 1" with infinitesimal generator C. The 
distribution of X(t) is given by v t . For s > 0, 

Ee~* z ? = J s- £ x < J] r(A,)M * 4 Sn (A)dA, 

i 

where the integral is along (upwards) vertical lines with !RA^ > for all i. 



The free energy for this model is given by [441 139] 

Urn ilog^ = inf[t -*(*)], 

n— >oo n *>0 

almost surely, where ^(z) — T' (z)/T(z). The conjectured KPZ scaling behaviour for the fluctuations of 
logZ™ was (essentially) established by Seppalainen and Valko |50j ; more recently, Borodin, Corwin and 
Ferrari [£Q [9] have proved the full KPZ universality conjecture for this model, namely that \ogZ™, suitably 
centered and rescaled, converges in law to the Tracy- Widom F2 distribution of random matrix theory. See 
also |51j . 

5. Reduced double Bruhat cells and their parameterisations 

The Weyl group associated with GL(n) is the symmetric group S n . Each element w £ S n has a 
representative w £ GL(n) defined as follows. Denote the standard generators for gl n by hi, ej and /j. For 
example, for n = 3, 

'0 0\ /0 0^ 

10, ^3=000 

v 0/ \0 I, 

o\ /o 0' 

h = 1 
0/ \o 

For adjacent transpositions Sj = (i, i + 1), define 

Si = exp(-ej) exp(/i) exp(-e 4 ) = (I - ei)(I + f t )(I - e,). 

In other words, Sj = ipi ^ where is the natural embedding of SL(2) into GLln) given by hi, 

and fi. For example, when n = 3, 
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Now let w = Sjj . . . Si r be a reduced decomposition and define w — . . . Si r . Note that uii = uv whenever 
l(uv) = l(u) + l(v). For n = 2, wo = s\ and 

For n = 3, wq — S1S2S1 = S2S1S2 is represented by 

u>o = = S2S1S2 = 

Denote the upper (respectively lower) triangular matrices in GL{n) by B and £?_, and the upper (re- 
spectively lower) uni-triangular matrices in GL(n) by N and 7V_. The group GL(n) has two Bruhat decom- 
positions 

GL(n)= (J B?LB = |J B-vB— 
ues„ ugs„ 
The double Bruhat cells G^ are defined, for u,v £ S n , by 

G"'" = BuBDB„vB_. 
The reduced double Bruhat cells L u '° are defined by 

I"- = NuN HB-vB-. 
We also define the opposite reduced double Bruhat cells M u ' v by 

M"'" = BuB n N^vN^. 



The reduced double Bruhat cell L w ' e (where e denotes the identity in S n ) admits the following parame- 
terisations, one for each reduced decomposition of w. See [361 [H l4l I16j. Set 

Yi(u) = (pi ( ™ „,-i), i = l, ...,n-l. 



I u 

Then, for any reduced decomposition i = %\ . . . i r of w, the map 

(ui, . . . ,u r ) n- Y h (ui) • • • Y lr (u r ) 

defines a bijection between C^ and L w,e . This bijection has the property that the totally positive part LQq 
of L u ' v corresponds precisely to the subset R> of C^ . There are explicit transition maps which relate the 
parameters (ui, . . . ,u r ) corresponding to different reduced decompositions of w. 

In the case n = 3, the two representations of an element in L Wo,e corresponding to the words 121 and 
212, denoting the corresponding parameters by (u±, 112,113) and (u^u^u^), respectively, are given by 

u lU3 \ (u' 2 

+ U2/U1U3 = \u' x u[u 3 /u 2 

1 l/u 3 1/M2/ V 1 u 3/ u 2 + 1 / u l l / u 'l u 3J 

The transition maps are given by 

(15) u[ = u 3 + u 2 /ui, u' 2 =uiu 3 , u 3 — u\u 2 /{u 2 + uiu 3 ). 

Their is a similar parameterisation for M w,e , due to Lusztig 1361 . For i = 1, . . . , n— 1, set Xi(v) — I+vfi. 
Take any reduced decomposition i = i\ . . . i r for w. Then the map 

(Wl, . . . ,V r ) H> Xi^Vi) ■ ■ ■ X ir (v r ) 

defines a bijection between C!^ and M w,e . This bijection also has the property that the totally positive part 
of M u,v corresponds precisely to the subset R> of CL. Q . 
In the case n = 3, the two representations of an element in M w °' e corresponding to the words 121 and 
212, denoting the corresponding parameters by (v\,v 2 ,v 3 ) and (v'i, v' 2 , v 3 ), respectively, are given by 

1 0\ / 1 

vi+v 3 1 = { v' 2 1 



v 2 v 3 v 2 1 1 \v[v' 2 v[+v' 3 1. 



with transition maps 



(16) U X = ; , V 2 =V!+V 3 , v 3 = ■ . 

Vl +V 3 Vi+ v 3 

We conclude this section with a simple lemma. Let b G G e ' w and write 6 = an where a = diag (eti, ...,a n ), 
say, and n G AT. Then [16j . for any w G S n , bw has a Gauss (or LDU) decomposition bw = [bw] _ [bw] [bw] + 
and nw has a Gauss decomposition nw — [nw] _ [nut] [nw] + . Moreover, [nw]_o = [niD]o G L w,e and 

G M™' 6 . Let i = i\ . . . i r be a reduced decomposition for u>. Then we can write 

M-0 = • ••iir( U r)) [M- = ^u( v l) • •••^trC'V)- 

Define Zj(u) = Lp i \ U ^-i)- Set a = a and, for 1 < k < r, a k = a k ~ 1 Z ik (u k ). Write a k — diag [a\, . . . , a*). 



^0 u 

Lemma 1. The following relations holds: [bw]o — a r and, for k = 1, . . . , r, Vk = uj~ l a^r^i/a^ 1 . 

Proof. Note that a[nw]_ = [bw]-o = [6u)]_[&u)]o, hence 

aY ix (ui) . ..Y ir (u r ) = X it (vi) . . . X ir (v r )[bw] . 

The result follows by repeated application of the identity aYi{u) = Xi(v)a', where a' — aZi(u) and v = 
u^ 1 a i+ i/a l . □ 



6. An evolution on upper triangular matrices 

As shown in [6], the path-transformations T w n can also be represented in terms of an evolution on the 
upper triangular matrices in GL(n,M). Let w — Si 1 ■ ■ ■ Si r be a reduced decomposition and n : (0,oo) — > K" 
a continuous path. Set tjq = rj and, for k < r, 

(17) n k = T ih . ..T h r) x k (t) = log / e- a ^"-^ds. 

Jo 

Then r\ r — T w r\ and, for each k < r, rjk — V + Ylj=i x j a ij ■ 

Write n(t) = (ijI, . . . , rf t L ). Define a path b{t) taking values in B by 

b i:j (t) =e 7),(t) / exp ( -^a fc (r?(s fc )) ] ds t ---ds^ 1 . 

iO<S J _i<S J _2 -<S ! <* \ k = i J 

If n is smooth, the b satisfies the ordinary differential equation 

(n n—1 \ 

53 h i dr l l + 53 e * dt 6 > 6 (°) = L 

i=l i=l / 

If n is a Brownian path (as in the next section) then b satisfies the equation interpreted as a Stranovich SDE. 
When n = 2, 

When n = 3, 

rf?? 1 dt 
<#>= I d*7 2 dt 



dry 
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and the solution is given by 



J _ J JQ<r<s<t 



b(t) = | e< f Q e<-<+rtds 

e''" 

Write b = an, where a = diag(e ??1 , . . . , e n ") and n G N. Set u k — e x * and Vk = e~~ Xk ~ ak< - Vk - 1 \ 

Theorem 3. [6j [7] -For eac/i t > 0, b(t)w has a Gauss decomposition bw = \bw\- [6iD]o[faZ'] + , wi£/i 
[6?D]o = exp(T to rj(/j)). Moreover, [nu)]_o = *ii ■ ■ • ii r («r) € L>q ■ 

By Lemma [TJ we also have = Xi^vi) ■ ■ ■ Xi r (v r ) S M™ Q e . 

6.1. The case n = 2. From the definitions: ai = e± — e2, wo = s i = s ei _ e2 , 



1 , 2 



a := ui =e Xl = / e^^+^ds w := v x = e^ 1 = e"" +r ' u^ 1 

ft 



o 



e 1 ? \ e 7 '' \ e 11 ' ) 1° 1 



'0 -V 

Taking wq = ( ^ ^ ) , we see that 



6w 

and 



1 OWe'u \ A -u- 1 
1 1 



nwo 



1 w\ (0 -1\ _ At -1\ _ At 0W1 -it" 
l) \1 ) ~ \l ) ~ \l u~ l ) \0 1 
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Hence 



[bw ] = e T ~o", [bw ]-=( 1 v fj=X 1 (v), [nteo]_o=(i u" 1 1 V| ' " ' 



as claimed. 

6.2. The case n = 3. From the definitions: 

oti = e-y — e 2 , a 2 = e 2 — e 3 , w = sis 2 si — s 2 sis 2 . 
For the reduced decomposition wq = s±s 2 sx, we have 



ui = e^ 1 = / ( 



u 2 = e X2 =l e-^+rfi^ds, =(e r,1 u 1 ,e r ' 2 u 2 /u 1 ,e ri3 /u2) 

us = e X3 
v x = e- yi = 



-1 


t 


e T)3 = 




















It 


1 e^-^+vl ds 


//o 


b = 





e" 2 






V° 









(e 1 ! 1 0\ /l 


Ml 




e" 2 


1 




v ef) \0 






1 



The identity 



follows from (fl5)) . Now, 



so we have 



e v 3 s -v 2 s +Vt d s = u 3 + u 2 /ui 



'0 V 

Wq = I -1 

.1 ; 



bwn = 



1 



e v u\u 3 



e r >'(u 3 + u 2 / Ul ) -e" 
e^ 3 

0\ (e^uius 



Ux + 1; 3 1 

«2 1/ 



e n u 2 /u\u 3 
\ 






e"7w 2/ 



'1 -I/U3 I/U1M3 



1 
.0 



-u 3 /u 2 - 1/ui 



and 



uiu 3 —ui 1 

1 0, 

U1U3 

"3 + "2/wi u 2 ju\u 3 



1 -1/U3 

1 



1 



l/«3 l/«2/ V° 



1/«1«3 

-w 3 /u 2 - 1/ui 



10 



Thus [bw ] Q = e T »o", 

1 0> 

[bw ]- = \ vi +v 3 1 ] = X 1 (v 1 )X2(v 2 )X 3 (v 3 ), 

V2V3 v 2 1, 

U1U3 

[m£j ]-o = ( "3 + u 2 /ui m 2 /wiu 3 ) = Yi(u 1 )Y 2 {u2)Y ? ,(u 3 ) 
1 l/u 3 l/it 2y 

as claimed. 

6.3. Evolution of the Lusztig parameters. As before, we introduce a probability measure V under 
which 77 is a Brownian motion in M" with a drift /1 and 77(0) = 0. For each k < n, set 

(r M ,...,T fe>fc ) = r( & V,...,r/). 

Note that this is given in terms of the principal minors b^ k \ k < n, of b by T^^rj 1 , . . . , f?' ) = log[M fe )u>Q fc ' ) ]o, 

(k) 

where Wq denotes the longest element in S k . The evolution of the triangular array Tk.j, 1 < j < k < n is 
given by (| 13|) . As remarked earlier, this process contains a number of projections which are also Markov. In 
particular, setting ^ k = T k _ k , we have, for k <n, 

d£ k = dr) k - e ik ' ik - 1 dt. 

This defines a simple interacting particle system on the real line which, in the coordinates £j and — 
1 < i < n — 1, has a product form invariant measure. There is an extension of this process, involving the 
Lusztig parameters, which is also Markov and, moreover, also has a product form invariant measure. Let 
v\, . . . , v q be the Lusztig parameters corresponding to a reduced decomposition wq — Si 1 . . . Si q , that is, 

[bw ]- = Xi^Vl) ■ --X^Vq). 

Set y k = — logWfc. The evolution of y k , 1 < k < q, is given by 

dy k = da ik (ri k ^i) + e~ Vk dt, 

where rj k = T ik . . .T n rj. Setting x k = y k - (*i k {ilk-i), note that dx k = e~ Vk dt and i] k = i] + J2j=i x i a j- 
Hence, 

k-l 

(18) dy k = da ik (r]) + ^ a lk {a i] )e~ Vi dt + e~ Vk dt. 

i=i 

Let /3i = on 1 and, for 2 < k < q, (3 k = Si 1 . . . Si k l oti k . Set 9 k — —(3 k {[i). If /i G wail = — f2, then 0^ > for 
all k and the diffusion has stationary distribution given by the product measure 

* = (g)r(6 k )- 1 ge k , 

fc=i 

where goidx) — exp(— 9x — e~ x )dx. This can be seen as a consequence of the following fact, which is the 
analogue in this setting of the output theorem for the M/M/l queue [44] . Let Xt be a standard one- 
dimensional Brownian motion with negative drift —9, and consider the one-dimensional diffusion 

dy = \/2dx + e~ y dt. 

This has a unique invariant distribution T(6)~ 1 gg. If we start this diffusion in equilibrium and define 
x t = Xt + 2(yo — yt), then i has the same law as x and, moreover, x s , s < t is independent of y Ul u > t, 
for all t. It follows that the measure tt is invariant. For an analytic proof of this fact, see |45j . See also [7j 
Proposition 5.9], where the equivalent property is proved in the 'zero-temperature' setting. 

If we choose the reduced decomposition i = 1 21 321 n — 1 n — 2 ... 21, and define, for m < n — 1 and 
1 < i < n — m, q m ,i = Ti +m> i + i - T 2+m _i^, then 



(2/1,2/2, ■ ■ • , y q ) = (51,1, 51,2, • ■ • , 9l,n, £2,1, ■ • ■ , 92,™-l, ■ ■ ■ , gn-i.i)- 

11 



Note that qn = ^+1 — £j, for 1 < i < n — 1. In these coordinates, the evolution is given by 



dq„ lt i — da t (i]) + e qm -'dt + YJ (2e~ 



,-<h,t+i 



e- q ^- 1 ) dt, 



1=1 



with the conventions that the empty sum is zero and q^o = +oo. Setting mi = n m +i — fi m , an invariant 
measure for this diffusion is given by the product measure (^) m i ge m 4 • The dynamics of this process can be 
viewed as a network, as follows. Consider the dynamics 

dQ = d(A -S) + e- Q dt, dD = dA- dQ, dT = dS + dQ. 

We think of A, S as the input and D 7 T as the output, and represent this system graphically as: 

A 

S 



D 



Then the evolution of the q m ,i can be represented as in Figure 2. To see directly from this picture the 
product-form invariant measure, note that, if A and S are independent standard one-dimensional Brownian 
motions with respective drifts A and a, with A < cr, then the diffusion Q has invariant distribution T(9)~ 1 gg, 
where 9 = a — A. Moreover, if we start this diffusion in equilibrium, then D t = A t + Qq — Qt and T t = 
St — Qo + Qt are independent standard one-dimensional Brownian motions with respective drifts A and a, 
and for each t > 0, (D S ,T S ), s < t is independent of Q u , u > t. The analogue of this fact in the setting 
of Poissonian queueing networks is the cornerstone of classical queueing theory. It is called the output, or 
Burke, theorem. Finally, we remark that the dynamics indicated by Figure 2 is the analogue, in this setting, 
of the dynamical interpretation given in |42j of the RSK correspondence as a kind of 'queueing network'. 



2,1 



T 3 ,i 





22,2 




?3,3 








91,2 




91,3 









?3,2 




92,1 




92,2 



93,1 



Ta;. 



Ta;. 



Taa 



Figure 2 . Graphical representation of the evolution of Lusztig parameters 



7. From the Feyman-Kac formula to Givental's integral formula 

The fact that the evolution equation (fT5|) for the Lusztig parameters has a product form invariant 
measure sheds some light on the relation between the Feyman-Kac formula (fit)]) and the integral formula of 
Givental. It follows from this that, for any given reduced decomposition of Wo, the random variables 

r ai ^ds, i = l,...,n-l 

can be expressed, via the transition maps, as rational functions of a collection of q = n(n — l)/2 independent 
Gamma-distributed random variables with respective parameters 9k, k < q, defined as above with ft = —rj. 
Note that f3 is a Brownian motion with drift A = — fj, G f2. Since the sets {9k, k < q} and {Xi — Xj, i < j} 
are the same, this allows (fTT))) to be written as a q-dimensional integral 

poo \ 

- E / e- Q '^)ds 

(19) =e A W / e-K le "° ,Wr '(" '.f]]^-^- '^, 

For example, when n = 3 and i = 121, we have 

0\ = Al — A2, #2 = Ai — A3, 03 = A2 — A3, 

and, using (fTT))) . 

, x 1 , 1 «i + u 3 

ri(Vl,W2,«3j = — , r 2 («l,U2,'U3j = — = ■ 

Vl V 1 V2V3 

In this case, the integral formula (TT9"|) becomes 

Aixi+A2a;2+A3a;3 / Ai — A 2 — 1 Ai— A 3 — 1 A 2 — A 3 — 1 



i* 



x exp [ — vi - v 2 — v 3 — e x i +;r2 J — e ^2+^3 1,1 + V3 j dvidv 2 dv 3 . 



Under the change of variables 



where T = (T^j, 1 < i < k < 3) is an array with (T 31 , T 32 , T33) = (xi, £2, X3), this integral becomes 

■0 A (a;) = / e ^i(' I 3i+T3 2 +T3 3 -T 21 -T 11 )+A 2 (T 21 +T 22 -T 11 )+A 3 Tii 

x exp (- e T32 - T21 - e T33 - T22 - e T22 - Tl1 - e T21 " T31 - e Tll_T22 - e T22 ~ T32 ) dT 11 dT 21 dT 22 . 

Since 'I' a (2;) is a symmetric function of A we see that this agrees with Givental's integral formula ([7]). 
We note that this is reminiscent of the derivation of Givental's formula given in |21j (see also |18L 119) )- 

8. Fundamental Whittaker functions 

The eigenvalue equation ([6]) also has series solutions known as fundamental Whittaker functions. Define 
a collection of analytic functions a nim (v), n > 2, m G (Z + ) n_1 , v G C" recursively by 

1 

a 2 , m {v) 



i\T(vi — V2 + m + 1) ' 
and for n > 2, 



n — 1 1 
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where [ii = i>i + v n j(n — 1), i < n — 1, and the sum is over k £ (Z+)" 2 satisfying fcj < m^, 1 < i < n — 2, 
with the convention that fco = = 0. Then |26[ Theorem 15] for each n, a n ^ m (y) satisfies the recursion 



n—X n — 2 n — 1 

III, 

i— 1 i— 1 i— 1 



i=l 



with the convention that a n , m = for to ^ (Z+)" , and a nj o(^) = IIi<j — Vj + 1) ■ Writing 
m ' = 2™=i m »( e i ~ e i+i), the series 

(20) TO,(:r) =]>>„,,„ Me-( m ' + ^ 

m 

is a fundamental Whittaker function as defined by Hashizume 25\ , and satisfies the eigenvalue equation ^ . 
We adopt a slightly different normalisation than the ones used in the papers [25] or [26] . Note that, for each 
x € M n , m v {x) is an analytic function of v. Moreover: 



Proposition 3. 



Proof. This comes from [3]. In the notation of that paper we are taking n = {cti/2, i = 1, . . . , n— 1}, 
m(2a) = 0, |?7 Q | 2 = 1 and ^ v {x) = 2<ik v {x) where q = n(n- l)/2. □ 

Now consider the function 6* t (x) defined by (fl4)) . Note that we can write 

g n (A) = — — i- — -fo(A) TT sin7r (^ ~ ^ _ 



Corollary 3. 



f27u)"n! 



9. Relativistic Toda and g-deformed Whittaker functions 

The algebraic structure underlying Theorem [T] is an intertwining relation between certain differential 
operators associated with the open quantum Toda chain with n particles. This structure should carry over 
to the setting of Ruijsenaars' relativistic Toda difference operators and g-deformed Whittaker functions \47\ 
\AS\ 115), 122] . A recent (related, but different) development along these lines is given in [8]. We will describe 
here the q- analogue of Theorem [T] in the rank one case, which corresponds to n = 2. 

In the case n = 2, the Whittaker function is given by 

i> x {x) - 2exp Q(Ai + \ 2 ){x x + x 2 f) K Xl _ X2 (2e^~ x 0/a) , 

where K u (z) is the Macdonald function. In this case, Theorem [T] is equivalent to the following theorem of 
Matsumoto and Yor [37] . 



Theorem 4. (1) Let (B^\t > 0) be a Brownian motion with drift fi, and define 



Then is a diffusion process with infinitesimal generator 

2^ + U l0g ^ (e ^Tx 

(2) The conditional law of B^ , given {Z^\ s < t;Z^ = z}, is given by the generalized inverse 
Gaussian distribution 

^(X I zY 1 e^ x exp (- cosh(a;)/z) dx. 
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Let < q < 1. Denote the g-Pochhammer symbol by (q) n — (q;q) n = (1 — q) • • • (1 — 9™) with the 
conventions that (q)o = 1 and (0)„ = 1. In what follows we also adopt the convention that 0° = 1. 
For A G C and z > 0, define 

» A(2„-s) 



This is a g-deformed Whittaker function associated with sh 22\. It satisfies the difference equation 

(i - q z+1 )Mz + 1) + M* - 1) = (g A + «- a )Va(*) 

where we set ip\{— 1) = 0, and is related to the g-Hermite polynomials by 



Fix 0<g<l,0<p<l and let (Y n , Z n ) n >Q be a Markov chain with state space {(y, z) E Z 2 : z > y > 0} 
and transition probabilities given by 

n((y, z), (y + l,z + I)) = p, U((y, z), (y, z + I)) = (I - p)q y , 

U(( y ,z),(y-l,z-l)) = (l-p)(l- q y). 

Note that V is itself a Markov chain with transition probabilities 

P(y,y + l)=p, P(y,y) = (l-p)qy, P(y,y - 1) = (1 -p)(l - 

and X = 2Y" - Z is a simple random walk on the integers which increases by one with probability p and 
decreases by one with probability 1 — p. Choose i/£M such that p = q v /(q v + q~ y )- 

Theorem 5. Let Yq = Zq = . The process (Z n , n > 0) is a Markov chain with transition probabilities 
Moreover, for each n > 0, the conditional distribution ofY n , given cr{Z m ,m < n} and Z„ = z, is given by 

K z (y) = i> v (z) — , y = 0,1,..., z. 

WyWz-y 

The proof is straightforward using the theory of Markov functions, by which it suffices to check that the 
transition operators II and Q satisfy the intertwining relation QK — KT\ where 

K(z,(y,z'))= 6 ^' q 



This intertwining relation is readily verified. When q = and v = 0, ^>u{z) — z and the above theorem can 
be interpreted as the discrete version of Pitman's '2M — X' theorem, which states that if X n is a simple 
symmetric random walk and M n = max m <„ X m , then 2M — X is a Markov chain with transition probabilities 
Q(z, z + 1) = (z + l)/2z, Q(z, z - 1) = (z - l)/2z. When q 1, it should rescale to Theorem|U 

The analogue of the output/Burke theorem in the setting of Theorem[5]is the following. If p < 1/2, then 
the Markov chain Y has a stationary distribution. If Yq is chosen according to this distribution and Zq = 0, 
the process (Z n , n > 0) is a simple random walk on the integers which increases by one with probability p 
and decreases by one with probability 1 — p. 
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